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We study chaos in a two dimensional Ising spin glass by finite temperature Monte Carlo simulations. 
We are able to detect chaos with respect to temperature changes as well as chaos with respect to 
changing the bonds, and find that the chaos exponents for these two cases are equal. Our value for 
the exponent appears to be consistent with that obtained in studies at zero temperature. 
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I. INTRODUCTION 

A characteristic feature of spin glasses is that the rel- 
ative orientations of the spins in the spin glass state are 
not uniquely determined by the model, but rather vary 
with external parameters such as the temperature or the 
magnetic field. At large separation, R, the correlation be- 
tween the spins varies in a chaotic manner as a function 
of temperature, and, when the temperature is altered by 
an amount AT, will change substantially at distances R 
greater than Zat where 
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which defines the chaos exponent £. This temperature 
induced chaos has been-.<gpite difficult to see in-Monte 
Carlo (MC) simulationarH or mean field theoryau, 
claims have been made that is absent or very smalfl 
A larger chaotic effect has been-pbserved when making a 
small change in the couplinggJ'B. 

Chaos with temperature has been shown to be a com- 
mon feature of the two niain models for the spin glas: 
phase: mean field theoryou and the droplet theoryErE 
The later is based upon real space renormalization group 
calculations which allow for quantitative results. In par- 
ticular, they indicate that a temperature perturbation 
generates a disorder perturbation and thus these twp per- 
turbations should have the same chaos exponentcE 

Here we study chaos with both AT and a change in 
the couplings, A J (defined precisely in Eq. (|J) below), 
by Monte Carlo simulations for a two-dimensional Ising 
spin glass at finite temperature. Since T c = for this 
model, all our data is in the paramagnetic phase, where 
the correlation function tends to zero at large distances, 
so we are looking for chaos in the sign of this decaying 
functioiD (such chaos has been shown in a one dimensional 
systems) . We study distances which are smaller than the 
correlation length so the chaps-we obtain is that corre- 
sponding to the critical pointMO. 

Our main results are as follows: 

1. It is possible to see chaos with AT as well as chaos 
with A J. 

2. The chaos exponents for AT and AJ appear to be 
equal. 



3. The chaos exponent found here at finite-T 
be consistent with that obtained at T = Ol 



ms to 



The plan of the paper is as follows: In §[jl]we define the 
model and various quantities of interest. § pTj| discusses 
finite-size effects which will be very important for the 
analysis while §IV presents the numerical results that 
are then interpreted in §|v|. 



II. THE MODEL 

We consider the Edwards- Anderson Hamiltonian with 
Ising spins and nearest neighbor couplings, 



H — JijSiSj 



(2) 



where the { Jy •} are drawn from a Gaussian distribution 
with zero mean and variance [J£] u equal to unity. We 
denote by [• ■ -] av an average over the interactions. The 
spins lie on a square lattice of linear size L with periodic 
boundary conditions. 

For each realization of the disorder we simulate several 
copies (or replicas) of the system. The basic quantity we 
calculate is the replica overlap 
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where a and b denote replicas and N = L 2 . When we 
investigate chaos with AT, some of the replicas will have 
identical bonds but slightly different temperatures and 
when we investigate chaos with A J some of the replicas 
will have slightly different couplings but the same tem- 
perature. 

Next we describe quantities that we calculated in the 
simulations. First of all, from replicas with the same 
temperatures and bonds, we compute the standard equi- 
librium quantities, g, the Binder ratio, and Xsg, the spin 
glass susceptibility, defined by 
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Xsg = L (q ) , 



(4) 
(5) 



where (■ • ■), denotes both the average over disorder and 
the statistical mechanics (Monte Carlo) average. During 
the simulation, the first to sweeps are used for equilibra- 
tion and the next to sweeps are used for measurements. 
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We check-that the system is in equilibrium by standard 
methods!!!}. The equilibration time to limits the maxi- 
mum size and minimum temperature we can study. In 
our case, we can reach T — 0.4 for L — 6 and T = 0.55 
for L = 10, which both require about to = 10 6 MC steps. 

Next we describe quantities that we calculate to deter- 
mine the chaos. First, consider chaos with changing the 
bonds, keeping the temperature fixed. This is done by 
running one replica with.a_set of bonds {Jij} and another 
with bonds {Jy}, wherelij 
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where Gaussian random variable with zero mean 

and unit variance. Note that the {J^}, and the {Jij} 
have the same distribution. A convenient measure of how 
much the change in the bonds alters the_spin orientations 
is the dimensionless "chaos parameter" B, 



(7) 



where the labels on the replicas refer to the bond distri- 
butions that are used. 

When the temperature is changed we consider the 
overlap from replicas at temperatures symmetrically dis- 
placed about T as follows: 
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where the temperatures are T± = T± AT/2. We believe 
that this is the first time that chaos with temperature 
has been calculated in this way. Other attemptsEl to look 
for chaos with temperature used a quantity which has 
an asymmetry in temperature and will therefore involve 
bigger corrections to scaling when the ratio AT/T is not 
vanishingly small. 

We average over from 80 to 400 realizations of the dis- 
order. Error bars are determined by grouping the results 
for the different samples into bias and calculating the 
standard deviations among binsli3. 
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FIG. 1. A scaling plot of the Binder ratio according to 
Eq. (^|) for different sizes and temperatures. The value of the 
correlation length exponent is v = 2.0. 



We also need the finite size scaling ansatzes for the 
chaos parameters, r AJ and r h 



r±j =f AJ (l Z AJ,L x / v t) 
r AT =f AT (l**AT, L l l v T\ , 



(11) 
(12) 



where we have used Eq. (Q) . It is inconvenient to analyze 
a function of two variables, so, following a suggestion of 
D. Huse, we have taken data where the second argument 
L~T, is roughly constant. We then try to collapse the 
data on to a single curve by plotting it against L^AJ 
with a suitable choice of £• 



IV. RESULTS 



III. FINITE SIZE EFFECTS 

To be in the scaling regime it is necessary to work 
at moderately low temperatures whereJinite size effects 
are important, and so finite-size scalingcj techniques are 
needed. 

Since £ ~ T _i/ , where £ is the bulk correlation length, 
and since, at the T = critical point, the ground state is 
unique, finite-size scaling predicts the following behavior 
for the Binder ratio and Xsg'- 

g = g {lVt) (9) 

Xsg = L 2 X SG (LV»T) . (10) 



Scaling plots for the Binder ratio and spin glass sus- 
ceptibility are shown in Figs. [I] and [|. In Fig. |l| the 
sizes L range from 4 to 12 and temperatures T from 0.4 
to 1.6. The resulting correlation length exponent at the 
T = transition is v = 2.0 ± 0.2. For Fig. | the sizes 
extend from L = 4 to L = 20, at temperatures from 0.4 
to 1.6, and the correlation length exponent is given by 
2/ =1.6 ±0.2. 

The value for v obtained from g is in good agree- 
ment jssiith finite temperature Monte Carlo simulations of 
Liaagtj and also agrees with the work of Kawashima et 
al.to. The estimate for v obtained from \ SG is somewhat 
smaller, presumably reflecting the systematic corrections 
to finite size scaling at the temperatures and sizes that 
we can study. Interestingly the same trend, namely a 
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FIG. 2. A scaling plot of xsg according to Eq. ( fj"o| ) for 
different sizes and temperatures. The value of the correlation 
length exponent is v = 1.6. 

larger value for v obtained, from g than from Xsg, has 
been seen in other modeldill. 

Combining the two values for v we estimate 



= 1.8 ±0.4 



(13) 



We also determined the spin-spin correlation function 
for temperatures between 0.8 and 1.2. By fitting this 
data to an exponential function of position, we estimate 
the correlation length, finding that it could be fitted to 
£ = AT~ V with A = 4 ± 0.5 and v = 1.8 ± 0.2, the latter 
being consistent with the estimates from the finite-size 
scaling analysis above. 

A scaling plot for chaos with AJ, following Eq. JTl| ) 
with L X I V T constant (and v = 1.8), is shown in Fig. [| for 
sizes between 4 and 10 with £ = 1.0. The perturbation, 
A J, lies in the range 0.05 — 0.3. Trying different values 
of £ we estimate 



C = 1.0 ±0.1 (chaos with A J) 



(14) 



A scaling plot for chaos with AT, following Eq. (|lz 
with L X I V T constant (and v — 1.8), is shown in Fig. 
for L < 10, and £ = 1.0. The perturbation, AT, lies 
in the range 0.05 — 0.4. Trying different values of £ we 
estimate 



C = 1.0 ±0.2 (chaos with AT) 



(15) 



Note that the exponents for chaos with AJ and AT, 
given in Eqs. (14) and (|l^), are equal within the uncer- 
tainties. We also see from Figs. |3| and ^ that the data for 
r AT does not deviate very much from unity, as compared 
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FIG. 3. A scaling plot of r A j according to Eq. ( |ll]) with 
[ J 1 / Vr p constant and v — 1.8, the average of our estimates 
from g and Xsg- The value of the chaos exponent is £ = 1.0. 
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FIG. 4. A scaling plot of r A T according to Eq. ( p^ ) with 
L 1 '"! constant and v = 1.8, the average of our estimates 
from g and Xsg- The value of the chaos exponent is £ — 1-0. 
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with the data for r AJ . This indicates that the amplitude 
of chaos with AT is smaller than that with AJ, even 
though the exponents are equal. 



V. DISCUSSION 

Studying the two-dimensional Ising spin glass by finite 
temperature Monte Carlo simulations, we are able to de- 
tect chaos with respect to both AJ and AT and show 
that the chaos exponents are equal, as expectedEI. 

We should point out that in order to define chaos with 
AT in the critical region it is necessary that v > To 
see this note that we need I at in Eq. (0) to be less that 
the correlation length, £ ~ T~ u , and also AT < T to be 
in the scaling region for chaos. In our case, this inequality 
is satisfied (note th at . ch aos in the critical region is also 
present when T c > 0E2fli3). Chaos with A J, on the other 
hand, can be defined irrespective of the relative values of 
v and 

Our estimates of the exponent, given in Eqs. (0) and 
(|l5|), are consistent with the value £ = 0.95±0.05Qfound 
from exact ground state determinations. A. similar value 
was also found earlier by Bray and Moored 

Finally-, we note that our value for v agrees with work 
of Lianglj who gets v ss 2 from Monte Carlo simulations. 
However a much larger value is inferred at T = Q_from 
domain wall renormalization group calculations^, i.e. 
v = 4.2 ± 0.5, and from exact ground state calculations, 
i.e. v = 3.559 ± 0.025. These discrepancies suggest a vi- 
olation of the scaling picture of the spin glass transition. 
Kawashima et al.EEl also find a discrepancies in the scal- 
ing theory. If there are violations of the scaling picture in 
two dimensions it would be very valuable to understand 
them since similar violations may also occur in higher 
dimensions with a finite T c , and also perhaps help re- 
solve disagreements between the droplet and mean field 
pictures. 



male Superieure, 69364 Lyon Cedex 7, France (Laboratoire 
associe au CNRS). 

1 F. Ritort, Phys. Rev. B 50, 6844 (1994). 

2 J. Kisker, L. Santen, M. Schreckenberg and H. Rieger, 
Phys. Rev. B 53, 6418 (1996). 

3 I. Kondor, J. Phys. A 22, L163 (1989). 

4 I. Kondor and A. Vegso, J. Phys. A. 26, L641 (1993); S. 
Franz and M. Ney-Nifle, J. Phys. A. 28, 2499 (1995). 

5 V. Azcoiti, E. Follana and F. Ritort, J. Phys. A. 28, 3863 

(1995) . 

6 H. Rieger, L. Santen, U. Blasum, M. Diehl and M. Junger, 
J. Phys. A. 29, 3939 (1996). 

7 A. J. Bray and M. A. Moore, Phys. Rev. Lett. 58, 57 

(1987) ; A. J. Bray, Comments Cond. Mat. Phys. 14, 21 

(1988) . 

8 D. S. Fisher and D. A. Huse, Phys. Rev. B 38, 386 (1988). 

9 M. Nine and H. J. Hilhorst, J. Phys. A. 24, 2397 (1991). 

10 M. Nine and H. J. Hilhorst, Phys. Rev. Lett. 68, 2992 
(1992); M. Ney-Nifle and H. J. Hilhorst, Physica A 193, 
48 (1993). 

11 One can also consider the opposite limit, R ^> £, in which 
case the chaos exponent ip-that for the infinite temperature 
paramagnetic fixed pointlij. 

12 R. N. Bhatt and A. P. Young, Phys. Rev. B 37, 5606 (1988); 

13 Also, as a check of the code, we verified that the Bethe- 
Peierls approximation, Xsg = (l + [w 2 ] a v) /(l - 3[u> 2 ] at >), 
where w = tanh(Jij/T), gives the same results as ours for 
T > 3. 

14 V. Privman, in Finite Size Scaling and Numerical Simula- 
tion of Statistical Systems, ed. V. Privman (World Scien- 
tific, Singapore, 1990), p. 1. 

15 S. Liang, Phys. Rev. Lett. 69, 2145 (1992). 

16 N. Kawashima, N. Hatano and M. Suzuki, J. Phys. A 25, 
4985 (1992). 

17 N. Kawashima and A. P. Young, Phys. Rev. B 53, R484 

(1996) . 

18 M.J. Thill and H.J. Hilhorst, J. Physique I 6, 67 (1996). 

19 D. Huse and I. Morgenstern, Phys. Rev. B 32, 3032 (1985). 



ACKNOWLEDGMENTS 



We thank Henk Hilhorst, Michel Gingras and David 
Huse for their comments on the first version of this pa- 
per. We appreciate Christophe Ney's help with C pro- 
gramming and thank Thierry Biben for sharing his com- 
puter. This research was supported by NSF Grant DMR 
94-11964. The work of MN was supported by a NATO 
fellowship and by the Centre National de la Recherche 
Scientifiquc. 



Permanent address: Laboratoire de Physique, Ecole Nor- 



4 



